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I. INTRODUCTION 



Several recent cosmological observations have revived the debate on the possibility of a 
cosmological constant A [|l]. The most discussed problem in this connection is the possible 
conflict between the age of the universe in the standard Friedmann- Robertson- Walker (FRW) 
cosmology, as inferred from recent measurements of the Hubble constant Hq 0, versus the 
age of the oldest stars in globular clusters . Constraints from the formation of large-scale 
structure in the universe together with dynamical estimates of its matter content, Q « 0.2 
Q, have also inspired suggestions of a non- vanishing cosmological constant, in order to 
make up fl + A /3Hq — 1 []5J. (Such a low value of fl may also be required to reconcile 
the observed high nucleonic content of clusters of galaxies, revealed through their X-ray 
emission, with the upper limit on the nucleon density from big bang nucleosynthesis (BBN) 
0.) These arguments limit the present value of the cosmological constant and the associated 
vacuum energy density, p v = A M|/87r, to be at most ~ (3 x 1CT 12 GeV) 4 , whereas quantum 
field theories and the inflationary paradigm generally demand a value for the vacuum energy 
before the start of the FRW cosmological evolution, which is larger by at least 100 orders 
of magnitude 0]. While there is, as yet, no theoretical understanding of this problem ||, a 
phenomenological approach is to invoke a time-varying cosmological "constant" . In such a 
scenario, the vacuum energy density decays from its original high value during inflation to 
the small value allowed for at the present epoch. Many such models have been proposed in 
recent years ||- |TB |. 

An important constraint on models in which the vacuum energy decays into radiation is 
provided by nucleosynthesis arguments, as was noted by Freese et al. [f| . Since their analysis, 
however, there have been several developments in this area [[17]]. Input parameters such as 
nuclear cross-sections, the neutron lifetime and the number of light neutrinos are now better 
known, while the theoretical calculation of elemental abundances has become more accurate. 
Moreover there have been major developments on the observational front, requiring revision 
of the abundance bounds adopted in earlier work. In this paper we correct and update 
this nucleosynthesis constraint and test the viability of phenomenological proposals for a 
decaying vacuum energy. 



II. BOUNDS FROM PRIMORDIAL NUCLEOSYNTHESIS 

In considering the possible form of the time variation of A(t), it is reasonable to assume 
that the vacuum energy is released as radiation, the energy density of which redshifts with 
the expansion scale factor as R~ A . (If the decaying vacuum energy creates massive particles 
instead, the present day energy density in matter will increase, exacerbating the cosmological 
age problem referred to above.) Hence Freese et al. |J considered that the vacuum energy 
density decreases proportionally as the energy density of radiation which they took to include 
e + e~ pairs and N v = 3, 4, 5 species of neutrinos in addition to photons, i.e. p T = p^+Pe+Pu = 
^g c sT A . To conveniently track the changing vacuum energy density they introduced the 
parameter 



2 



which they supposed remains constant with the expansion, even as p v decays increasing p r . 
Although apparently plausible, this assumption is in fact physically inconsistent during the 
nucleosynthesis era for the following reasons. Since the weak interactions freeze out at a 
few MeV, the neutrinos generated by the decaying vacuum energy cannot be thermalized at 
lower temperatures. Further, e + e~ pairs turn non-relativistic and annihilate soon thereafter. 
Consequently, to have p v decrease proportionally to T 4 as desired for phenomenological 
reasons, we must assume that the vacuum energy decays into photons alone and redefine 

x = . (2) 

P 7 + Pv 

This does remain constant throughout the BBN epoch, since the photons produced by the 
decaying vacuum energy do get thermalized (within a Hubble time) down to a temperature 
of about 10 keV. 

To quantitatively study the effects on nucleosynthesis, two important changes have to 
be made to the BBN computer code [13]. First, we add a term to the Friedmann equation 



to account for the vacuum energy density: 



57T 



(We have neglected the energy density of non-relativistic matter, p m , as is appropriate in 
the radiation-dominated era.) Secondly, the equation of local energy conservation now has 
a term p v describing the production of energy by the decaying vacuum energy: 

R 

Pv + P + 4-p = . (4) 

Freese et al. [0] had stated that this leads to the modified time-temperature relationship 

T(£)oc[s cff (l-x)r 1/4 t- 1/2 ; (5) 

this, however, obtains only if the effects of e + e~ annihilation as well as the relative decrease 
of the neutrino temperature due to vacuum decay are neglected. (Because of the entropy 
generated by vacuum energy decay, the temperature of the background neutrinos (as well 
as the nucleon-to-photon ratio if) continually decreases in ratio to the photon temperature 
even at T-Cm e and drops below its standard asymptotic value of T u = (4/ll) 1 / 3 T.) We 
discuss this in the Appendix and implement these effects in the Wagoner code as updated 



by Kawano [18]. We assume 3 light neutrino species, use the latest value of the neutron 



lifetime r n = 887 ±2 s |JL9[|, and incorporate small corrections to the helium abundance as 



reviewed in ref. [[L7|] 



The overall effect of a vacuum energy decaying to produce entropy is a decrease in the 
synthesized abundance of 4 He. This is because the effect of the increased expansion rate, 
which causes earlier neutron-proton freeze-out, is sub-dominant relative to the effect of the 
dilution of rj, which delays the onset of nuclear reactions. Consequently, although the n/p 
ratio at freeze-out is higher, it decreases by /3-decay to a smaller value before the neutrons 
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can be synthesized in light nuclei.^ Thus the magnitude of the decaying vacuum energy 
during nucleosynthesis can be bounded by requiring that the abundances of the synthesized 
elements be within observational limits. The relevant bounds considered in earlier work |9| 
were 

0.22 < F p ( 4 He) < 0.26, 10~ 5 < (D/H) p < 10" 4 , (6) 

for the primordial mass fraction of helium and the primordial abundance by number of 
deuterium, respectively. We adopt these bounds as well in order to illustrate the difference 
in results due to our use of an improved code and compare the constraints obtained on x as 
a function of r\ in Figure |TJ(a), having rescaled the previous result || to correspond to our 
redefinition fl2|) of x, taking p r = yP 7 , for N v = 3. (Note that the constraints corresponding 
to the upper bound on Y p are out of range, to the left of the y-axis). The constraint quoted 
earlier |{J corresponds to a bound on the vacuum energy which is less restrictive than our 
result by a factor of about 2. 

Next we impose the updated bounds on the element abundances inferred from recent 



observations of helium in metal-poor extragalactic HII regions pC§ and of deuterium in high 
redshift clouds along the line of sight to quasars f2T|j: 



0.23 < F P ( 4 He) < 0.25, 1.1 x 10~ 5 < (D/H) p < 2.5 x 10" 4 . (7) 

(We do not consider the bound on 7 Li because it does not provide a useful constraint.) We 
emphasize that these are conservative bounds based upon consideration of a variety of data 
which are critically discussed in ref. |]T7[. Because of the uncertainties in the input nuclear 
cross-sections and the neutron lifetime, as well as residual numerical errors in the computer 
code, the yield of an element for some specified value of rj and x is not uniquely determined 
but has a spread which is about ±0.5% for helium and as much as ±50% for deuterium. 
Using a Monte Carlo method in which the BBN code is run many times while the input 
parameters are sampled at random from their known distributions, the "95%" CL. bound 
can be identified as the value at which 5% of the runs pass the imposed observational bounds 



22| . On the basis of these results we find a maximum value of 

x max = 0.13 (8) 

for i] ~ 3.7 x 10~ 10 as shown in Figure |l|(b). This corresponds to the bound 

p v < 4.5 x 10" 12 GeV 4 , (9) 

evaluated at a temperature T = 2.6 MeV. This is about 3 times more stringent than the 
often quoted constraint calculated by Freese et al. |J. For comparison, a constant vacuum 
energy is constrained by nucleosynthesis to be p v < 4.4 x 10 -17 GeV 4 for n in the range 
10~ n - 10~ 9 . (Although the bound on a time-varying cosmological "constant" appears 
weaker, it should be noted that it would decay to a much smaller value by the present epoch 
than one which is truly constant.) 



x Note that if the decays do create (electron) neutrinos as assumed by Freese et al. Q, then there 
will be further direct effects on neutron-proton interconversions; such effects, however, were not 
considered by these authors. 
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FIGURES 





x 

FIG. 1. Bounds on the decaying vacuum energy parameter x = p v / (/> 7 + p v ) as a function of the 
nucleon-to-photon ratio i] (evaluated at T = 10 8 K). The shaded regions are allowed. In panel (a) 
we compare previous results (dashed lines) with those obtained using the improved nucleosynthesis 
code, adopting the same bounds on the elemental abundances. In panel (b) we adopt updated 
observational bounds and also show (dashed lines) the "95%" C.L. limits corresponding to the 
known uncertainties in input parameters. 
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III. CRITIQUE OF DECAYING A MODELS 



Many of the proposed models in the literature assume that 

A oc R~ m , (10) 



where 1 < m < 3 |]10| . Since the associated vacuum energy density p v then decreases more 
slowly with time than the energy density of radiation, even a p v having its maximum al- 
lowed value today would have been negligible at the epoch of nucleosynthesis. However, 
it would be physically more consistent to normalize the vacuum energy density at some 
high energy scale instead of at the present epoch. For example Gasperini [|ll|] interprets the 



cosmological constant as a measure of the Hawking temperature of the De Sitter vacuum, 
T v = (A/127T 2 ) 1 / 2 ; he assumes that the vacuum temperature equals the radiation tempera- 
ture close to the Planck era, but that the subsequent expansion causes T v to decay faster 
than T so as to be consistent with the upper bound on the cosmological constant today. 
However, the value of p v during nucleosynthesis would still be substantial; in order to obey 
the constraint ([8]), the exponent in eq.(^) is now required to be 

m>3.5. (11) 



Thus Gasperini's suggestion is inconsistent with many models [|TUJ which assume m to be 
less than 3. (Note that it does not matter that the constraint (||) was derived assuming 
m = 4 since it would be even more stringent for smaller m.) 



Some other models [|T^| have too many free parameters and/or are not specified suffi- 
ciently explicitly to be confronted with BBN. Specific models motivated by physical con- 
siderations fl3|-|T6l| can however be constrained, either by the bound (||), or by the related 
bounds on the expansion rate or the values of fundamental constants during BBN. In some 
cases the authors have already considered such constraints but commented that the associ- 
ated uncertainties preclude a definitive test. We reexamine the question taking into account 



all such sources of error and imposing conservative constraints from BBN ||17|| . 



Lima and collaborators Jl3| consider nonsingular deflationary cosmology models with a 



decaying vacuum energy density. The cosmological history starts with the decay of a De 
Sitter vacuum and evolves smoothly to a quasi-FRW stage at late times. Irrespective of the 
spatial curvature, the models are characterized by the time scale iff 1 which determines the 
initial temperature of the universe as well as the largest value of the vacuum energy density. 
The present matter and radiation content of the universe is generated by the slow decay of 
the vacuum energy density; for its time dependence, the authors use the phenomenological 
ansatz 

am p (, 1-PH\ 

Pv = _E =M ^ 1 + _£_j , ( 12 ) 

where p to t = p v + Pm + Pv and (3 is a dimensionless parameter of order unity. This auto- 
matically generates a primordial inflationary scenario at the time iff . The model thus 
represents a generalization of ref. || since it does not fix the spatial curvature to be zero 
and introduces a time dependence in the parameter 
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x' = - P - =p + (l-p)¥ r . (13) 

(Pv+Pm + Pr) Hj 

At late times, when the term H / H\ is negligible, we recover the model by Freese et al. 
and find that x' ~ (3. In the physically preferred case, the deflationary process begins at the 
Planck time, Hf^Mp , while the latest epoch at which it can occur is electroweak symme- 
try breaking (to permit the subsequent generation of a baryon asymmetry), i.e. iff^M^ 1 , 
so H/H\ will necessarily be negligible during nucleosynthesis. Thus the nucleosynthesis 
constraint (^) translates into the requirement 

(3 < 0.13 ; (14) 

the inequality is reinforced by the fact that x as defined by us (eq.(||)) is always bigger 
than x'. This conflicts with the lower limit on (3 [|13j following from the requirement of a 
sufficiently long age for the universe: 



(3 > 0.21 , (15) 



thus excluding this model IT 



The model proposed by Wetterich |TJ| involves a scalar field with the potential 




Focexp [—£-04 , (16) 



where a > is a parameter. Such exponential potentials may well arise in models of unifica- 
tion with gravity such as Kaluza-Klein theories, supergravity theories or string theories. The 
scalar field couples to gravity and its energy-momentum tensor induces a time- dependent 
cosmological "constant" . The model further contains a coupling of the scalar field to matter 
cx (3p, where (3 is another model parameter which may, in principle, have different values 
during the matter-dominated and radiation-dominated eras. Wetterich assumes that (3 X is 
zero and finds that (3 m < from considerations of the stability of condensed objects with 
a static energy density. After some critical transition time t tr , the cosmological "constant" 
adjusts itself dynamically to become proportional to the energy density p (of radiation plus 
matter), and all energy densities (p, V, 2 /2) decay oc t~ 2 . Wetterich considers two possible 
scenarios for the initial epoch before this transition. In the first, the universe is initially 
p-dominated and the scalar field contributions do not influence the time evolution of the 
scale factor. In the second, the universe starts off being ^-dominated; this would lead to 
an exponential increase of the scale factor as in vacuum energy driven inflation, although 
without any subsequent entropy production.^ As an alternative to the above possibilities, 
a purely scalar-gravity system is also considered; this necessitates additional fine-tuning to 
evade the severe experimental bounds on such a theory. Thus only the first case, viz. initial 
p- domination, appears to be well motivated. 



2 This is unappealing since the horizon and flatness problems of the universe cannot be solved if 
there is no such entropy production p3| . 
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The parameter (3 can now be constrained as follows. In the matter-dominated era for 
distances much smaller than H -1 , the coupled system of gravity and small scalar fluctuations 
around the cosmological value <f>(t) is identical to the standard Brans-Dicke (BD) theory. 
Experimental bounds on the Brans-Dicke parameter u then constrain /3 m to be small 

\P m \ < 0.016 . 

If in fact (3 m is identically zero, then the contribution of the scalar field today is 



(17) 



n 



<t>0 



3 

2^ 



i-a 



(18) 



Thus flfo > corresponds to a 2 > 3/2. We can obtain a stronger bound on a if the transition 
time ttr is smaller than the BBN epoch so that the increased total energy density during 
BBN due to the scalar contributions speeds up the expansion as t — > t' = £ _1 t. The recently 
revised [[HJ upper limit on such a speed-up factor £ is 



£ < 1.12 



(19) 



corresponding to the conservative abundance bounds in eq.([?|).f] This translates into the 
lower bound 



a 2 > 9.9 



where we have used the relation £ 2 — 1 = 2/|a 2 — 2|. In turn this requires 



n mo > 0.85 



(20) 



(21) 



at the present epoch, i.e. the present cosmological constant contribution must be quite 
small. A possible escape route is to assume f3 m ^ 0, but this leads to a time dependence of 
the nucleon mass and, possibly, other fundamental constants (which can also be constrained 
by nucleosynthesis if the model is made explicit enough). Even so the modification to the 
age of the universe, to — \Hq 1 {1 — £), is negligibly small, removing the major motivation 
for the scenario. This conclusion may however be evaded by assuming that the parameter a 
is not constant but varies with the curvature scalar or the value of the scalar field. Another 
escape route would be to assume that t tr is larger than the BBN epoch, so the expansion 
rate at this time is unaltered. It is evident that the requirement of successful nucleosynthesis 
imposes interesting constraints on the possibilities considered in ref. (TJj . 
Moffat 



15 



considers a model in the context of a BD theory where the non-minimal 
coupling of a scalar field to gravity allows for a negative pressure associated with the kinetic 
energy of the field. The cosmological "constant" is assumed to be space-time dependent and 
is equated to the kinetic energy in the BD theory. As a result the evolutionary equation 
for the scalar field in a FRW universe produces an attractor mechanism which drives A 
towards a minimum of the potential and leads to a small constant value for A today. In this 
model Moffat finds a solution to the problem of the age of the universe with Q v several times 



3 Wetterich [14| takes this limit to be £ 2 — 1 < 0.1, which then requires a 2 > 22. 
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larger than Q m at present. However, the constraint (|19D on the speed-up factor £ during 
nucleosynthesis requires Q v ^ 0.1, thus removing the motivation for this model. 

Berman [|16[] assumes A oc t~ 2 in a BD theory, leading to a possible time variation in the 
gravitational constant Gn- If, in fact, Gn is constant, the vacuum energy density equals the 
energy density in radiation, so that x as defined in eq.(|l|) is 0.5 and thus in conflict with 
the nucleosynthesis constraint On the other hand, if Gn varies with time, it must do 
so in this model as Gn oc t^ 1 during matter-domination and as Gn oc t 2 during radiation- 
domination |1| . Such a strong time variation is ruled out given that the value of Gn during 



nucleosynthesis is required to be within about 25% of its present value 



IV. CONCLUSIONS 



We have examined proposed models of a time-varying cosmological "constant" which can 
significantly increase the expansion age of the universe, even for a large Hubble parameter 
today, thus evading the potential age problem. However, these models also imply a deviation 
from the standard expansion history during primordial nucleosynthesis. We find that even 
conservative observational limits on the abundances of the synthesized light elements imply 
severe constraints on such models, demonstrating the power of this cosmological probe. 
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APPENDIX: 



As mentioned earlier, the ratio of neutrino to photon temperature continues to decrease 
even after e + e~ annihilation. To study this, we follow the standard procedure |25| of com- 
bining the equations describing the change in entropy and energy conservation to obtain 



S(T) - S(Ti) 



4ir 2 x 



15 Z - 1 J Ti 



dT R 3 T 2 



(Al) 



where p v = xp 7 /(l — x) and p 7 = 7r 2 T 4 /15. Inserting the expressions for S(T) and S(Ti), 
S(T) = R?{p + p)/T, we can rewrite the above equation to obtain 



rp3 



1 + 



45 
2^ 



r21 



03 



+ 



12x 



11(1 -x)R 3 Jt z 



df R 3 T 2 , 



(A2) 



m e /T 



dy y r 



m e 



n/2 



1 



ey + 1 



(A3) 



The first term in the sum above tracks the standard behaviour of the neutrino temperature 
during e + e~ annihilation. The numerical solution of this equation is shown in Figure § for 
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several values of the parameter x; note that the curves converge at the neutrino decoupling 
temperature. 

Freese et al. || did not take into account the differing neutrino and photon temperatures 
or the change in the number of degrees of freedom due to e + e~ annihilation, hence their 
relation for the nucleon-to-photon ratio, rj = 251 oc T Zx ^ l ~ x \ only holds after this epoch. 
Between neutrino decoupling and e + e~ annihilation, the correct relation should be 

7] oc T l2x/ll{1 ~ x) . (A4) 

Similarly, the equation for R(T) as given by Freese et al, viz. R oc T 1 /^'- 1 ), which holds 
only at T<^m e , is modified to 

R oc rCH-T^/ii^-i) (A5) 
between neutrino decoupling and e + e~ annihilation. 
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